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Abstract 

We describe effects of anisotropy caused by the crystal lattice in d-wave 
superconductors using effective free energy approach in which only one or- 
der parameter, the d-wave order parameter field, is used. All the effects of 
rotational symmetry breaking, including that of the s-wave mixing, can be 
parametrized by a single four derivative term. We find solutions for single 
vortex and the vortex lattice. Extending the formalism to include the time 
dependence, effects of anisotropy on moving vortex structure are calculated. 
Both direct and Hall I-V curves as functions of the angle between the current 

and the crystal lattice orientation are obtained. 
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It is widely believed that the superconductivity in layered high Tc cuprates is largely 
due to the d(^x2_y2^ pairing |T|. It has been observed that this influences the vortex lattice 
structure [0,0. There are also indications that even though the major bulk pairing mecha- 
nism is of d wave nature, there is a small admixture of the s-wave pairs in the condensate. 
Ren et al 0, using a phenomenological microscopic model (in the weak-coupling limit), and 
Soininen et al 0, considering attractive nearest neighbors interaction, proposed an effec- 
tive Ginzburg-Landau (GL) type theory using two fields. Two order parameters, s and d, 
describe the gap functions in corresponding channels. The free energy was constructed to 
include the fourfold symmetry 0- 

+7,|ns|2 + 7d|m|2 ^ ^^[5*(nj _ n^)^ + c.c] (i) 

where 11 = —iV — e*A (with c = h = 1). Using equations following from this free energy 
or using more fundamental equations §], one obtains a characteristic four-lobe structure for 
the s-wave inside a single vortex ||7|J9| jTI|] . The distribution of the magnetic field was also 
obtained recently j^J^. The s-wave vanishes outside the core, while the d-wave distribution 
becomes rotationally invariant, indistinguishable from the usual Abrikosov solution. Since 
the fourfold vortex core structure is seemingly in conflict with the symmetry of the triangular 
lattice, the asymmetry of vortices therefore can, in principle, distort the usual triangular 
vortex lattice. If one were to look for differences in the behavior of vortices between this case 
and that of the conventional s-wave superconductors, one would like to be closer to Hc2, so 
that the influence of the core will be more significant. Another phenomenon in which the 
core plays a major role is the dissipation in the course of the flux flow. 

In this paper we study the above two phenomena using time independent and time de- 
pendent effective GL equations within a simplified one component framework. The simplicity 
of the formulation allows us to clarify, essentially without loss of generality and the use of 
numerical methods, various delicate questions about single vortex and the vortex lattice. 
The degrees of freedom we include in the analysis contain: (1) an arbitrary rotation angle (p 
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between the crystal lattice and the vortex lattice and (2) all the possible lattices, not only 
the rectangular ones considered before PJ^,p!l[1). This is the first time that the lattice is 
demonstrated to be centered rectangular ( CR) with the most general lattice included in the 
analysis. Moreover the treatment can be extended to moving flux lattices, which, as is well 
known |T2|JT^ are more demanding, as far as calculational complexity is concerned. 



Within the two-field formulation, Soininen at al observed that, in a predominantly 
d-wave superconductor, the s-wave component is generally very small: it is "induced" by 
the variations of the larger d component. In the bulk, only the d field acquires a nonzero 
value, while, near the core, the rotationally noninvariant gradient term s*(n^ — n^)(i "com- 
municates" the deviations from the condensate value of d to that of s. Even near the core, 
where d is small, it is still larger then s by a factor of 20. Since the s field is "induced" by d, 



it varies on the scale of the d-wave coherence length = y7d/«d- Then 7sn^s ~ {.Is/Qt) ^ 
is small compared to agS if (75/7^) (ad/as) ^ 1 (typically 7^/7^ ~ 1). This is an excellent 
approximation in both near and far regions from the core |]7|,p!0|,p| . Therefore, the field s is. 



to the first order in l/os, s = — (7t,/as)(n^ — Ill)d. Substituting this back to GL equations, 
we obtain, to first order in l/a^, the effective free energy 

^'"^"^^ = i'""^'' ~ + - vd* (nj - Ulfd. (2) 

Here we have replaced by a more conventional notation l/2md and defined 77 = j^/as- 
The last term should be treated perturbation. 

Note that the two-field equations are highly nontrivial even at the linearized level, so 
that authors of Ref. resort to the variational estimate in order to find a solution. The lin- 
earized one-field equation can be solved perturbatively in t]. The advantage of this equation, 
especially as far as relation to experiments is concerned, is that the number of coefficients is 
much smaller: instead of 10 parameters in the two-field free energy, there is just 3. One can 
further motivate the use of the effective free energy Eq.(|^) even with no connection to the 
two-field formalism. Generally superconductivity is a phenomenon of "spontaneous gauge 
symmetry breaking " phenomenon irrespective of the mechanism of pairing or channels in 



which it occurs. The symmetry should be represented by a single order parameter: the su- 
perconducting phase. While in pure s-wave or d-wave superconductors the phase is simply 
identified with the phase of the gap function, in more complicated microscopic theories with 
few channels opened, the superconducting phase is just the common phase of various gap 
functions. Therefore, quantities other than the phase which enter various phenomenologi- 
cal GL type equations, although useful, are not directly related to the spontaneous gauge 
symmetry breaking. 

The general form of the effective free energy can be obtained just by the dimensional 
analysis and symmetry. It is well known that in a symmetric field theory, at the level of 
the dimension three (relevant) terms, full rotational symmetry is restored. Therefore, one 
has to consider "irrelevant" dimension five terms in order to break the rotational symmetry 
down to -D4/1. d*{Ily—Il'^)'^d is the only such term which does not respect rotational symmetry 
and therefore is important for studying anisotropy effects. Consequently only one coefficient 
7] is needed to describe the anisotropy effects. The origin of rj might, in principle, come from 
sources other than the d-s mixing. The one field formulation also avoids the problem of the 
artificial second phase transition at Tg assuming = ot'{Ts — T) in the two-field formulation. 

Single vortex solution near Hd- In strongly type II materials, we can safely ignore 
the magnetic field. The GL equation can be solved perturbatively as follows: Write d = 
do + Xdi , where do = /o(?")e*'^ is the solution of the standard unperturbed GL equation and 
A = irjm'^ad is a dimensionless small parameter. The angular dependence of di is easily 
found to contain only three harmonics e"^*"^, e'^'^'^ and e^**^ 

rfi(r,0) = f_,(r)e-'^^ + Mr)e'^ + f,ir)e''t (3) 

This is consistent with the fourfold symmetry which is built into the theory. The analytic 
expression for /q does not exist, but there are a number of known approximations. Using 
one of them, fo{r) = r / ^Jr^~+~^ , the set of linear equations for fi, f2, fs are then solved 
numerically (the /s equation decouples from the other two) ||TB[. The d-wave configuration 
is nearly indistinguishable from that of the two-field formalism as expected, while the s-wave 
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configuration calculated with s = — (7„/as)(n^ — Ill)d has the four-lobe structure with four 
zeroes on the axes (same as in 0] but different from |10[| ). 

Vortex lattice. In the opposite limit near Hc2 one first neglects the nonlinear terms in 
the GL equation and finds the lowest energy solutions "^kni^^v) of the linearized equation. 
The lattice solution is constructed from the linear superposition of \E^fc^(x,?/) in such a 
way that it is invariant under the corresponding symmetry group of the lattice: d{x, y) = 
J^nCn'^knix^y). A general lattice in 2D can be specified by three parameters a,b and a, 
where a and b are the two lattice constants, while a is the angle between the two primitive 
lattice vectors as defined in Ref. [|14|. Flux quantization provides a constraint: Hab sin a = 
$0- In the d-wave superconductors the rotational symmetry is broken, therefore the relative 
orientation of the vortex lattice to the underlying lattice becomes important. We denote (p 
to be the angle between a and one of the axes of the underlying lattice. The most general 



lattice is built into the solution via appropriate relations among C„ |TJ] . The lattice structure 
is then found by minimizing the free energy or equivalently, the Abrikosov's parameter 



/3a 



Here (...) is the average over space. 



We solve the linearized GL equation perturbatively in the anisotropy parameter rj in the 
Landau gauge, A = {0, Hx){see details in [|r§]): 



ipix,y) 



1 



1/4 

2 ■ exp{iyk„ 
H / 



1 + 7]mde*H—Hi — 
16 \Ih 



X exp 



-^-(^-kl'^' 



(4) 
(5) 



where H4 is the fourth Hermit polynomial and Ih = yl/e*H. The upper critical field is: 

2m(ia' r 



(T, - T) + 8r]mja' (T, - T) 



(6) 



Note that Hc2{T) bends upwards around Tc, in agreement with the two-field results [0. 
This effect has been reported in some experiments. However one should be cautioned 
against taking this too seriously. As we have discussed, there is another rotationally invari- 
ant term of the same dimensionality r(i*(H^)^(i which gives contribution similar to Eq.(^): 
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AHc2{T) = — r ^"^^" (Tc — T)^ . This means that for large enough positive r the sign of the 
-f^c2 (7") curvature changes. 

To present Pa it is convenient to use variables: p + ia = ( = ^exp(i«) ||14|. After some 
rather nontrivial calculations we obtain, = + r]'(3\ with 



cr 



E MO 



+ 



E Wn+1/2(C) 



(7) 



/3\ = ^V^ReU'^ 



E<'(C) 



5:^„(C)G„(a) 



+ \ n n + -,n 



^+2 



where r/' = rjmde*!! is the dimensionless parameter appropriate for this situation, m„(C) = 
exp(27riCn^) and G'n(o") = {QAn'^a^n'^ — ASTran'^ + 3). Having calculated the Abrikosov pa- 
rameter Pa, one finds the vortex structure by minimizing it with respect to ip, p, and a. The 
minimization with respect to the angle (f is easily done analytically, while further minimiza- 
tion of /3™"(p, cr) is done numerically. The angle a as function of rj' is given on Fig.l. The 
corresponding if is zero. For a = 74° one gets rj' ~ 0.02. There is a phase transition from CR 
to square lattice at rj'^ = .0235. Note that this calculation, unlike that for the single vortex, 
is valid for arbitrary n. Using standard methods, one can take into account variations of the 
magnetic field and calculate corrections to the magnetization curve using the Pa calculated 
here [0 . 

Moving lattice. The time dependent GL equation describing the time evolution of the 
order parameter the one-field formulation is 



7 + te*^) d = - (^n^ -ad)d + r]{Ul - UD^d - 2p\d\% 



(8) 



where $ is the electric potential. It involves just one additional parameter 7 compared to the 
2x2 matrix for the two-field formalism. Although, in principle, this parameter, describing 
various dissipation effects, can have a complex part we will consider only real values. 
To generalize the above procedure for finding the structure for a moving vortex lattice near 
Hc2, one considers the motion caused by an electric field E making an angle with the 
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crystal [1,0,0] axis. The vortex lattice velocity is: E = — v x B. For a general direction of 
the electric field the fourfold symmetry of the system is completely (explicitly) broken. 

Even for the simple s-wave the problem of finding the moving lattice solution is nontrivial. 
However in that case there exists the "Galilean boost" trick |13| to solve the linearized 



problem, which is not applicable to the d-wave. Our approach is to use perturbation theory 
in rj to find a complete set of solutions of the linearized equations and then impose the 
periodicity conditions to construct the vortex lattice. It is more convenient to perform 
the first step in the gauge aligned in the direction of the electric field, in order to make 
both the scalar and vector potentials independent of y and t. For the second step however, 
it is preferable to use a gauge aligned in the direction of the vortex lattice. However in 
general the vortex lattice will not be periodic along this special direction. To construct 
this general periodic solution, one has to solve a periodicity constraint equation for the 
coefficients C^, where k is now a continuous index. We then combined the two steps using 



gauge transformation [0. The results are as follows. The upper critical field of course 
depends on velocity 

iic2iT) =ho + h,{T, -T) + h,{T, - Tf (9) 

with ho = (m^e*) [-1 + (9 + cos46)r7m|)7^f^] 7^t;^, hi = 2a'{md/e*) {1 - 12r]mlj'^v'^), 
h2 = IGa'^rjm^/e* . Note that the curvature /12 hasn't changed compared to the static case, 
but we have two new effects. First of all, the electric field (or, equivalently, electric current) 
shifts Hc2 by a negative constant (proportional to E"^ ) to a lower value. It simply depends on 
the angle . This is expected. Secondly, although the curvature h2 doesn't change compared 
with the static case, the slope hi acquires a negative contribution proportional to E"^. In 
the s - wave case the phase boundary was first found and discussed in Ref. [0]. There are 
a couple of peculiarities associated with it like the existence of a metastable normal state 
and the unstable superconductive state. The same applies to the present case. As far as 
we know, these peculiarities haven't been directly observed in low Tc materials. It would be 
interesting to reconsider this question for the high Tc materials. 



After a lengthy and nontrivial calculation |T5[, the vortex lattice solution is amazingly 
simple: 



1 fH\'/^ 



exp ( iknv) exp 



X 



1 + 



3im{6+(/p) 



with Ci 



1 + e 



m=l 



C2 



-f^m 1 ^n^H - ige 



-^(l + 3e 



(10) 



-4j0 



^ige-^''', C4 



\^g-4ie g^j^j^ ^ ^ rjj^g standard Abrikosov's procedure to develop an approximation 
for small order parameter around ifc2 can be applied also in the flux flow case (see [0]). 
Using this expression the correction term in the expansion of the Abrikosov parameter (3a 
changes and now the function G'„(cr)in Eq.(|^) depends in a very simple way on electric field: 
G„((t) = e^*^(647rV2n^-487ran2 + 3)-8e2*'^^2cos2e(87r(Tn2-l). One immediately observes 
a surprising fact - the dependence on the angle G and velocity g is only via the combination 
g^ cos 20. For example, the resulting lattice for O = 7r/4 and arbitrary g will be the same as 
without electric field at all! Also apparent complete breaking of the rotational symmetry by 
general direction of the electric field is not felt by 13a- Indeed, the lattice for some arbitrary 
O and g"^ is the same as for = and g^' = cos 20. However, the fourfold symmetry 
has been reduced. For fixed rj' and (7^ cos 20, the minimization was performed numerically 
and we again obtain only CR lattices. The angle a turns out only weakly depend on the 
combination (7^ cos 20 . 

This is the first time that the lattice is demonstrated to be CR type using the most 
general lattice in the analysis. It turns out that general oblique lattices have higher energy 
than the CR ones despite the fact that rotational symmetry is completely broken by both the 
electric field and by the underlying atomic lattice. Intuitively in the symmetric case this fact 
is understandable: CR lattices are more symmetric, however for rotationally nonsymmetric 
or moving lattices this is no longer so. 

I - V curves for the flux flow. Now we consider the dissipation in vortex cores due to flux 
flow. As it is well known, the fourfold symmetry forces the conductivity tensor cTjj, defined 
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by Ji = CijEj, to be rotationally symmetric. Namely, aij = aSij + a^Sij. Here a is the usual 
(Ohmic) conductivity, is the Hall conductivity and Sij is the antisymmetric tensor. The 
additional term in the free energy corrects the values of a and , but the correction is of 
the order rj and therefore small. So, to see anisotropy, we definitely would like to go beyond 
linear response. This has been done for simple s - wave TDGL [0| near ■ We will neglect 



pinning and consider motion of a very large bundle. While there is a normal component of 
the conductivity, here we will concentrate on the contribution of the supercurrent only. For 
the discussion of the relative contribution of the two see [l^ . 



We therefore calculate the current as a function of E beyond linear response. Note that 
in this case the condensate has to be properly normalized: < d*d >= a^/ {2PPa)- Using the 
normalized d, we found the anisotropic direct and Hall currents 

AJ^.. = ^^^(l + cos4e), (11) 
A3HaU = -V^-0^sin4e. (12) 

Note that both direct and Hall currents depend only on the fourth harmonics of the angle 
between E and the crystal lattice orientation. Only the cubic power of E contributes, all 
the higher orders terms are cancelled. 

To summarize, we studied the static and dynamical anisotropy effects in d-wave super- 
conductors using the effective one component model in which these effects are parametrized 
by a single parameter 7]. For static lattice we clarified several issues and found the critical 
value of 7]'^ = .0235 at which thermodynamic transition to the square lattice occurs. (The 
sample of [0] is found to be in CR phase while that of [Q] is in square phase). The results 
of our study of the moving lattices for the direct current. Hall current magnetization and 
phase boundary are very simple and are given in Eqs.|rT],|T^,|^ and ^ respectively. It would 
be interesting to determine rj' of the same sample from a few of the above effects. 
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Hu, C.S. Ting and F.C. Zhang for discussions and correspondence. The work is supported 
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Figure Caption: 

Fig. 1 The angle a of the centered rectangular lattice as a function of the parameter 77' 
describing the asymmetry. There is a phase transition from centered rectangular to square 
lattice at 77^ = 0.0235. 
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